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Abstract 

Using the Madelung fluid picture of Schrodinger equation for a single 
free particle moving in one spatial dimension, we shall specify a class of 
wave functions whose quantum probability density is being trapped by the 
quantum potential it itself generates. The global convexity of the quantum 
potential generated by the initial self-trapped wave function will then be 
shown to further localize the quantum probability density through a self- 
generated focusing equation for a finite interval of time. 
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1. Motivation: Schrodinger Equation as Self- Referential Dynamics 

Let us consider a quantum system of a single free particle with mass m 
moving in one dimensional space, q. The wave function ip{q\t), where t is 
time, is then governed by the linear Schrodinger equation 

ihdtij{q-t) = -^dli,{q-t). (1) 

Quantum theory claims that the quantum probability density, p(g) = 1^(0' 
gives the essential information about the position of the particle [H, 0, 0, 
Hence, it is important to study the dynamics of p{q;t). 

One way to keep track directly the dynamics of p(g; t) is to use Madelung 
fluid hydrodynamical interpretation of the Schrodinger equation To do 
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this, let us write the wave function in polar form, ip{q; t) = R{q; t) exp{iS{q; t)/h), 
where the quantum amplitude R and the quantum phase S are real-valued 
functions. Inserting this into Eq. ([1]) and separating the real and imaginary 
parts, one gets 

a.^ + %^ + t/ = o, 

2m 

d,p + d,[^p)=0, (2) 

where U is the so-called quantum potential generated by the quantum prob- 
ability density as 

Next, let us define a velocity field as 

v{q-t) = ^d,S{q-t). (4) 

Taking the spatial gradient to both sides of the upper equation in ([2]), and 
using Eq. (jlj) one has 

m% ^ -a,U. (5) 

Due to its similarity with Euler equation, the term on the right hand side of 
the above equation, —dgll, is usually called as quantum force field. Moreover, 
again using Eq. (jlj), the lower equation in ([2]) can be seen as the continuity 
equation for p{q;t). 

One of the great lesson from the Madelung fluid picture is its self-referential 
property. Namely, the quantum probability density p{q; t) will generate a 
quantum potential U {q; t) which in turn will rule the quantum probability 
density the way it must flow and so on and so forth. There is a dynamics 
circularity between the object to be ruled and the rule itself. Self-referential 
property is a simpton of complex non-linear systems, and is argued to be the 
general origin for the emergence of many interesting phenomena [6] . We can 
then expect to find interesting self-organized behavior difficult to be observed 
directly from Eq. ([T]), at least in a certain time scale. 

To see an example of the effect of the self-referential property of the 
Madelung fluid and to help describing the motivation of the present work. 
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let us discuss the dynamics of a Gaussian wave packet. In this case the 
quantum probabihty density will take the form [3] 

PiT,t) = le.v(-^), + (6) 



Here, Z is the trivial normalization constant, at is the width of the wave 
packet at time t and a is its width at t = 0. The corresponding quantum 
potential can thus be calculated to give 

U{q;t) = -^^+-^. (7) 

One can clearly see from the global concavity of the quantum potential 
given above that the quantum force field is repelling the quantum probability 
density to spread out. As time goes, the quantum probability density keeps 
its Gaussian form but is rescaled to be wider. This will also rescale the 
quantum potential to be wider while keeping its form, thus keeps the effect 
of repelling the quantum probability density. This simple picture is missing 
in the original Schrodinger equation of ([1]). Hence, Gaussian wave packet 
provides an example of a class of wave functions being repelled by its own 
self-generated quantum potential. 

The question that we want to discuss in this paper is then: is there any 
class of wave functions whose quantum probability density is self-trapped or 
self- attracted by the quantum potential it itself generates? 



2. Self- Trapped quantum probability density 

Obviously, the self-trapped quantum probability density that we want to 
develop must decay as the self-generated quantum potential increases. One 
therefore has to develop another relation between p and U which provides 
the supplement to Eq. ([3]). There are many possibilities to create an p — U 
relation which satisfies the above mentioned condition. In this paper, we shall 
consider a class of wave functions {^}, whose quantum probability density 
is given as the exponential of its quantum potential as follows 

P(g;i) = ^exp(-/?[/(g;t)), (8) 

where /3 is a positive real- valued contant, and Z is a normalization factor 
given hj Z = f dqexp{—pU). It is then clear that p decays as U grows. 
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Figure 1: The profile of quantum potential U{q) (solid blue line) satisfying Eq. ([9]), and 
the corresponding quantum probability density p{q) (dashed line) which generates the 
quantum potential. The solid red curve is the approximation of U{q) given by Eq. (|lip . 

Notice that Eq. ([8]) together with the definition of quantum potential given 
in Eq. Q comprise a differential equation for p(g) or U{q), subjected to the 
condition that p{q) must be normalized. 

In term of U, taking the spatial derivative twice to both sides of Eq. ([H]) 
and recalling the definition of quantum potential given in Eq. one obtains 

d',U = ^{d,Ur + A% (9) 

where A = ^^/Am/h?f3. First, one can see from the above equation that 
the positivity of quantum potential, U > 0, guarantees its local convexity, 
dqll > 0. Further, it is clear that the above differential equation is invariant 
under the transformation q —q. These facts will guide us later to find a 
class of solutions which are positive everywhere (thus convex everywhere), 
and symmetric with respect to the vertical line g = 0. 

Let us first solve Eq. ([9]) for a spatial region very close to the extremum 
points of the quantum probability density or quantum potential, Iq' — Q| -C 1, 
where Q is the extremum point under consideration satisfying dqU{Q) = 
dqU\Q = 0. In this case, assuming the continuity of dgU, the first term on 
the right hand side can be ignored, dqU ~ 0. Hence, in this spatial region 
we have 

dlU^A% |g-g|<l. (10) 

The above differential equation can then be solved with constraint dqUlq = 
and U{Q) = Uq, to give 

f/(g) ^f/QCOsh(±A(g-g)), |g-Q|<l. (11) 
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The point Q turns out to be the local minimum (maximum) point of the 
quantum potential U{q) if Uq is positive (negative). 

To check the desirable trapping property, one must see the global behavior 
of the quantum probability density and its corresponding quantum potential. 
To do this one should solve equation Eq. ([9]) for the whole space. We shall 
resort to numerical methods to do this. For the reason mentioned previously, 
we shall confine ourselves to a class of solutions which satisfy the following 
boundary conditions dqU\Q = and Uq > 0. Namely, we want to find a class 
of solutions where the first spatial derivative at the origin is vanishing, and let 
Uq to varies only to positive values. The first condition will explore the sym- 
metry of the differential equation, while the second boundary condition will 
guarantee that the quantum potential is nowhere non-convex. Note that the 
case f/o = is trivial. Namely, imposing the boundary condition dqU\o = 
to Eq. ([9]) one gets d'^U\o = such that U{q) = everywhere. Hence, p{q) 
is unnormalizable. This case therefore is excluded from our consideration. 
Fig. [1] shows the numerical solution of Eq. iQ with the boundary condition 
dqU\o = and Uq = 1. The profile of quantum probability density (dashed 
line) is plotted together with the profile of quantum potential it itself gener- 
ates (solid blue line). For comparation, we have also plotted the approximate 
solution given in Eq. (ITT]) (solid red line). 

First, it is observed that very near to the extremum point, q = 0, Eq. (ITT!) 
fits the numerical solution very well. Yet it immediately departs from the 
numerical solution as |g| is getting larger. We can see clearly from the numer- 
ical solution that globally the quantum probability density is indeed being 
trapped by the quantum potential it itself generates. This global property 
can be understood for any positive value of f/o > as follows. First from Eq. 
(jnD, inserting the boundary conditions at g = 0, one gets dgU\o = MJq > 0, 
such that U{q) is locally convex at g = 0. Hence, since dgU\o = 0, at spatial 
points nearby g = one has U{q) > Uq > 0. Moreover since the second term 
on the right hand side of Eq. (Q is always non-negative, at this region one 
has d'^U > 0. This geometrical reasoning can be extended such that for the 
whole space one gets 

d'^U > 0. (12) 

Namely, the case when Uq > will give an everywhere positive and convex 
quantum potential. Hence, q = turns out to be the global minimum of 
U{q). Using this fact in Eq. ([8]) one can conclude that quantum probability 
density, p(g), is being trapped by its own self-generated quantum potential. 
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U{q). 

Next, from the definite positivity of the quantum potential and Eq. ([3]), 
one has 

(9ji? < 0. (13) 

Hence, the quantum amphtude is everywhere concave. Supphed with the 
fact that R{q) is finite, one can conclude that R{q) must cross the g— axis 
at finite value of g = ±qm- This means that the support of R{q), [—qmyQm] 
is finite. At g = ig^, f/(g) is infinite, yet the quantum probability density 
is vanishing. Further, notice that at g = ig^, dqR{q) is discontinuous. 
However, dqp{q) = 2R{q)dqR{q) = is continuous. Hence, the self-trapped 
wave functions is a truly localized wave function with finite tail (support). 

Now let us compare the class of self-trapped wave functions we developed 
in this paper with Gaussian wave packet. As discussed in section [H the latter 
represents a class of self- repelled wave functions. Fig. [2] shows a self-trapped 
quantum probability density (solid line) plotted together with the Gaussian 
wave packet (dashed line) of the same second moment. One can see that 
Gaussian is more localized in the region near the maximum. Yet, in contrast 
to the self-trapped quantum probability density which possesses support only 
on a finite region, the Gaussian wave packet possesses support on the whole 
space. 

Another interesting difference between Gaussian and self-trapped wave 
function is that though both are quite similar to each other, their corre- 
sponding quantum potentials are very contrast from each other. While the 
quantum potential generated by Gaussian wave packet is concave everywhere, 
the quantum potential of self-trapped wave function developed in this paper 
is convex everywhere. As will be discussed in the next section, one can then 
expect that this makes their future destiny very different from each other at 
least for a finite period of time. 

3. Self-Generated Focusing Equation 

Intuitively, for spatially one dimensional case we are considering, the 
global convexity of the quantum potential makes the quantum force field 
to drag the quantum probability density toward the center, i.e. the global 
minimum of the quantum potential, at least for some interval of time. To 
see this more clearly, let us take the spatial derivative to both sides of Eq. 
([5]) to give us 

mdt9 = -m0^ - ^Jf/. (14) 
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Figure 2: The profile of a self-trapped quantum probability density (solid line) compared 
with a Gaussian quantum probability density (dashed line) with the same second moment. 

Here dt = d/dt, 9 = dqV = dgS/m is the velocity divergence and we have 
used the following fact from the differential calculus dt = dt + vdg. Hence, 
for the class of self-trapped wave functions we are considering here, where 
d^U > everywhere, one will get 

dtO < 0, (15) 

which is a focusing equation. Let us note that the above relation is valid 
for any value of 0{q), hence for any form of quantum phase S{q). One can 
thus conclude that in contrast to Gaussian wave packet which is self-repelled 
everywhere, the self-trapped wave function is self-focused everywhere. 

Due to the continuity of the evolution of U {q; t) which is guaranteed by 
the continuity equation in the lower part of ([2]), U{q;t) will keep its global 
convexity for a certain finite interval of time, say T. Given an initial self- 
trapped quantum probability density, p(g;0), T of course depends on the 
initial quantum phase as well. Thus it depends on 6^0, T = T{9o), where 6*0 is 
the value of 9 = {l/m)d'^S at t = 0. Yet, the finite- ness of T is independent 
of the choice of 9q. Namely, no initial wave function can avoid the focusing 
equation of ( |T5l) . In particular, during the interval of time T one has 

dt9 + 9^ = -{l/m)d'^U <0. (16) 

This will imply dt{9^^) > 1, such that one obtains 

9'\t)>9^^+t. (17) 

The above inequality shows that if 9q is negative, then within a time < 
l/\9o\, 9-^ must pass zero, namely, one has 9 —oo. Of course this will 
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occur as long as 

l/^o<T(^o), (18) 

posseses a solution. Hence, caustics will develop in a finite interval of time. 
During this time, the quantum probability density keeps converging toward 
the origin. 

4. Conclusion and Discussion 

To conclude, first we have specified a class of wave functions whose quan- 
tum probabihty densities are being trapped by the quantum potentials they 
themselves generate. The self-trapped wave functions are then dragged by a 
self-generated focusing equation of f|T5l) to localize even further, at least for 
certain interval of time depending on the initial quantum phase. 

For a given initial self-trapped quantum probability density, if there is Oq 
such that inequality ffTSl) is satisfied, then the velocity divergence is equal to 
minus infinity everywhere, 6 — * — oo, in a finite time. During this period of 
time, the quantum probability density converges everywhere until caustics is 
developed. The next question is then to investigate whether the singularity 
in velocity field (caustic) will also induce a singularity in quantum potential 
and how the quantum probability density looks like at this time. It is also 
tempting to study what happens to the wave function after the development 
of caustics. We leave this for future work. 
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